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Abstract. In this communication we report on a peculiar property of barrier 
transmission that systems governed by the nonlinear Schrodinger equation share with 
the linear one: For unit transmission the potential can be divided at an arbitrary 
point into two sub-potentials, a left and a right on(\ which have exactly the same 
transmission. This is a rare case of an exact property of a nonlinear wave function which 
will be of interest, e.g., for studies of coherent transport of Bose-Einstein condensates 
through mesoscopic waveguides. 
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Quantum mechanics is full of surprises, in some cases even in quite elementary 
situations, that are seemingly well known from textbooks. A recent example is the 
following remarkable observation for barrier penetration. 

Let us consider the transmission of a one-dimensional wavefunction through a 
potential V{x) with V{x) — >■ for x — » ±00. A well known example is a rectangular 
potential well 



(Vo > 0) with a transmission probability [1] 

(m is the mass of the particle and E the energy). At certain "resonance" energies Ej-^s 
the potential is 100% transparent, for example at En = —Vq + (/i^vr^/Sma^) ra^ > 0, 
n = n^ain, nmin + 1, • • •, for the rectangular well ([T]). Such resonances have been 
studied for more general potentials, as for example symmetrical or asymmetrical double- 
barrier structures in connection with resonant-tunneling, where they are denoted as unit 
resonances (see, e.g., [2] and references therein). 

Chabanov and Zakhariev [3] discovered that, at such a resonant energy, the 
potential scattering shows a surprising symmetry: Dividing V{x) into two distinct parts 
at some arbitrary point x' defining a 'left' and a 'right' potential 

J ^(^) a; < x' \ X < x' , . 

^^(") = | x>x' ' ^^(") = | V(x) x>x' ' 

the transmission probabilities |TlP for Vi, and |TrP for Vr are equal at the resonance 
energy, 

|TLp(Eres) = |TR,|2(Eres), (4) 

despite the fact that the left and right potentials can be very different. This is shown 
easily using the transfer matrix M connecting the amplitudes of the wavefunction on 
the left hand side, ipi^) = Aexp{ikx) +Bexp{—ikx), with those on the right hand side, 
ip{x) = C exp(i/cx) + D exp{—ikx) (the limit — > oo is understood for a potential not 
vanishing outside a finite range): 

(5) 

where the unit determinant detM = |ap — = 1 guarantees that the reflection 
probability |-Rp = and the transmission probability |Tp = sum to unity 

(see, e.g., [1]). If the potential is cut at x' according to (jsj), the total transfer matrix 
can be expressed as a product of the transfer matrices Mr and Ml of the right and left 
potentials, respectively: M = MrMl. For unit transmission we have |a| = 1 and 

P = ttR/^L + /5r«l = (6) 
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which imphes /^R/aR = — /3£/aL and hence equal reflection probabilities of Va and Vl. 
This property, which can also be extended to more general situations, offers an intuitive 
approach to the design of potentials with desired transmission properties (see [3] and 
references therein for more details). 

In this communication we report an even more puzzling fact, namely that these 
properties remain valid for the nonlinear Schrodinger equation (NLSE) 

|!_^" + (^_V)^_^|^|2^ = 0. (7) 

The nonlinearity destroys the superposition principle in linear quantum mechanics and 
we therefore cannot base our proof on matrix techniques as above. The NLSE ([T]), also 
known as Gross-Pitaevskii equation, attracted much interest in recent years because it 
describes the dynamics of Bose-Einstein condensates in a mean field approximation at 
low temperature (see, e.g., [4]). Note that the chemical potential /i takes over the role 
of the energy E. We furthermore note that also the nonlinearity g may be x-dependent. 

Transport properties of cold atomic gases in designed mesoscopic waveguides are 
of recent interest [5] . Here in particular barrier transmission of the nonlinear waves has 
been studied in a number of papers [5-8] that assume an experimental setup in which 
matter waves from a large reservoir of condensed atoms at chemical potential fi are 
injected into a one-dimensional waveguide in which the condensate can propagate. In 
these articles it was shown that the results obtained from the stationary NLSE ([T]) are 
in excellent agreement with numerical solutions of the time-dependent NLSE 

ihij{x,t) = -—i:"{x,t)+V{x)ij{x,t)+g\ij{x,t)\^ij{x,t) + foexp{-i^t/h)6{x-Xo) (8) 
2m 

where the source term foexp{—ifit/h)6{x — Xq) located at x = Xq emits monochromatic 
matter waves at chemical potential n and thus simulates the coupling to a reservoir. 
The barrier potential V{x) is assumed to be zero for x < Xq. As discussed in these 
articles, the definition of a transmission probability is not unambiguous, because of 
the implicit dependence on the scattering wave function |'?/'(a;)p. Different prescriptions 
assuming either a fixed source strength /o or a fixed incoming current can be formulated 
which differ in a strongly nonlinear situation. No ambiguity appears, of course, if the 
nonlinearity vanishes in the region far away from the scattering potential. Note that 
for unit transmission the differences between fixed output or input conditions as well as 
the definition ambiguities disappear. See [5, 8] for more details. 
A fixed output boundary condition 

ip{x) Cexp{ikx), ip'{x) ikC exp{ikx) for x — > +oo (9) 

with k = ^J2m{^ — g\C\'^)/h determines the wavefunction il){x) and its derivative il)'{x) 
in the region — oo < x < oo uniquely [9]. This provides a recipe to define an effective 
incoming amplitude A in the upstream region and a transmission coefficient 

|T|2 = (A;/A;a)|C/A|2 with k^ = ^2m{fi - g\A\^)/h . (10) 

For the fixed output problem, the transmission coefficient is a unique function of the 
chemical potential /i. If the problem is formulated as a fixed input problem, one observes 
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a bending over of the transmission spectra and an interaction induced bistability for 
strong nonlinearity where the transmission coefficient is no longer a unique function of 
/i. In addition, it was found that the points of unit transmission, the resonances /ires, 
survive in the nonhnear case g 0, of course shifted to different positions depending 
on g. In the following examples, we will use transmission coefficients derived from the 
fixed output condition 

If the potential V{x) is 100% transparent at the chemical potential fi = /^res, the 
solution in the far upstream region is given by 

ip{x) C exp(ifcx + lip), ip'{x) — >• ikC exp(i/ca; + 'up) for x — > — oo . (11) 

where (y9 is a real- valued phase. 

Now we define a second fixed output problem given by the NLSE 

+ {^l-V)x-g\xW = ^ (12) 

with the initial conditions 

xi.^) C* exp(— ifca; — iip), x'{^) ^ —ikC* exp{—ikx — iip) (13) 

for X — > — oo, which corresponds to scattering through V{x) from the opposite side of 
the well. Noting that, for real values of /x, the complex conjugate of a solution also 
solves the NLSE, we find immediately that the solution of (12) and (13) is given by 

X{x) = V'*(x), x'ix) = ip'*{x), (14) 

and in particular 

x{x) = C* exp(— i/cx), x{^) = —ikC* exp(— ifcx) for x +oo . (15) 

This means that the barrier potential V{x) is also transparent for waves with chemical 
potential fi coming in from the opposite side of the barrier even if the potential is 
asymmetrical. Note that this is not satisfied automatically, because the equality of the 
left-to-right and right-to-left transmission probabilities valid in the linear case is not 
guaranteed for the NLSE, at least to the knowledge of the authors. 

Now we divide the barrier potential V{x) as given in equation ^ and obtain two 
new fixed output problems: The first one is determined by the NLSE 

+ (/i - ^r)^r - gli^nUn = (16) 

with the initial conditions 

iPr{x) C exp(iA;x), ip'ji^x) — » ikC exp(iA;x) for x +oo , (17) 

which corresponds to tunneling through the potential V^^l^x) from the left to the right 
hand side of the barrier. The second one is determined by the NLSE 

XL + (/"-^L)XL-^7|XLrXL = (18) 



2m 

with the initial conditions 



Xl{x) — » C* exp(— i/cx — iip), x'l ~^ —ikC* exp(— iA;x — 'up) (19) 
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for X — oo, which corresponds to tunnehng through the potential from the right 

to the left hand side of the barrier. 



From equation ( 14 ) we obtain at the cutting position x' 

Xl{x') = ipU^'), Xhix') = -ipRix') . (20) 

Since the wavefunctions in the respective potential-free regions are uniquely determined 
by the values of the wavefunction and its derivative at x' it follows together with the 
initial conditions (17) and (19) that the asymptotic wavefunction Xl{x) for x ±oo 
is equal to the asymptotic wavefunction ip^ii^) ^ -Foo up to an irrelevant phase 
factor exp{iip). Thus the respective transmission coefficients coincide no matter which 
particular definition of the transmission coefficient is used. 

For illustrative purposes we explicitly demonstrate in the following the equality 
of the two transmission coefficients for the case that a fixed source strength /o is 
assumed in equation ([s]). In section II. A of [8] it was shown that the source strength 
is connected with the effective incoming amplitude via /o = I^/crAr with = 
A/2m(yU — By considering the stationary solutions of equation (8| it was 

further shown that A-^ is determined by the wavefunction and its derivative at the 
position xo of the source. Choosing xq = x' we obtain from equation (12) in [8] 

2ifcRAR = tp'^ix') + ik'M^') (21) 

with k' = A/2m(/i — g\iljR{x')\'^) / h = ^y2m{jj, — glxhix')]"^) / h. Since Xl{x) corresponds 
to scattering from the opposite direction the respective equation for the effective 
incoming amplitude differs by a sign. With the definition /cl = A/2m(/i — g\Ai^\'^)/h 
we thus arrive at 

2ifcLAL = x'dx') - iA;'xL(x') = iij'^ix') + ik'M^')T = (2iA;RAR)* (22) 

using (20) and (21). Consequently, we obtain = |^rP and /cl = so that the 

respective transmission coefficients 

\Tl\' = = = \Tk\' (23) 



coincide. More precisely, the left-to-right transmission of Vl is equal to the right-to-left 
transmission of Vr (compare the remark above). 

As an illustration, we present results based on a numerical solution of the NLSE for 
two cases. Figure [l] shows the transmission probability for a rectangular well ([T| with 
m = h = 1, Vq = 50, a = 20 and an attractive nonlinearity g = —1 acting along the 
whole X-axis. The outgoing boundary conditions are chosen as C = 1. We observe unit 
transmission, |Tp = 1, at the resonance /ires ~ 2.135. 

We now divide the rectangular well with width 2a = 40 at the value x' = —10 into 
two separate parts and calculate the transmission coefficient of each of these rectangular 
wells Vl (width 10) and (width 30) separately. Their transmission curves displayed 
in figure [T] cross at the chemical potential /^res, exactly as predicted. 

As a second example we choose the double-Gaussian barrier 

V{x) = Vo [exp(-(x + bf/a^) + exp(-(x - bf/a^)] (24) 
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Figure 1. Transmission coefficient (-.-.) of the NLSE for the rectangular potential 
well ([ij (parameters m — h — 1, depth Vq — 50, width 2a — 40 and attractive 
nonlinearity g — —1). Also shown are the transmission probabilities for the left (width 

10, ) and right potential (width 30, — ) obtained from dividing V{x) into two 

parts. These two probabilities are equal at the resonance at /ires where the potential 
V{x) is transparent. 



with parameters Vq = 1, b = 7.35, a = 6/5 and a nonlinearity of g = 0.005 studied 
in [6,8], where the apparently small value of g nevertheless causes a strong deviation 
of the nonlinear resonance curve from the linear one. In figure [2] we observe full 
transparency at /i ~ 0.7632, i.e. below the potential maxima. In addition the figure 
shows the transmission curves for the left and right potentials, where the potential is 
divided hj x = —8.0. Again the left and right transmission coefficients are equal at the 
point of transparency. 




Figure 2. Same as figure [T] however for the double barrier potential (24) 



In conclusion, we have shown that the interesting property of linear quantum barrier 
penetration at unit transmission, the equality of transmission through the arbitrarily 
divided sub-potentials, can be carried over to the nonlinear case which is much less 
understood. In addition of being of interest in its own right, this remarkable fact may 
provide an intuitive way to treat transport in nonlinear quantum systems, as for example 
the Bose-Einstein condensates in the mean-field approximation. 
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